Constraints on Cold Dark Matter Accelerating Cosmologies and Cluster Formation 
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We discuss the properties of homogeneous and isotropic flat cosmologies in which the present ac- 
celerating stage is powered only by the gravitationally induced creation of cold dark matter (CCDM) 
particles {Q.m ~ !)• For some matter creation rates proposed in the literature, we show that the 
main cosmological functions such as the scale factor of the universe, the Hubble expansion rate, 
the growth factor and the cluster formation rate are analytically defined. The best CCDM scenario 
has only one free parameter and our joint analysis involving BAO + CMB + SNe la data yields 
Sim ~ 0.28 ± 0.01 (la) where f2m is the observed matter density parameter. In particular, this 
implies that the model has no dark energy but the part of the matter that is effectively clustering is 
in good agreement with the latest determinations from large scale structure. The growth of pertur- 
bation and the formation of galaxy clusters in such scenarios are also investigated. Despite the fact 
that both scenarios may share the same Hubble expansion, we find that matter creation cosmologies 
predict stronger small scale dynamics which implies a faster growth rate of perturbations with re- 
spect to the usual ACDM cosmology. Such results point to the possibility of a crucial observational 
test confronting CCDM with ACDM scenarios trough a more detailed analysis involving CMB, weak 
lensing, as well as the large scale structure. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X 
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1. INTRODUCTION 

The analysis of high quality cosmological data (e.g. su- 
pernovae type la, CMB, galaxy clustering, etc) have con- 
verged towards a cosmic expansion history that involves 
a spatially flat geometry and some sort of dark energy 
in order to explain the recent accelerating expansion of 
the universe [11-0I- The simplest dark energy candidate 
corresponds to a cosmological constant, A (see (sl-fioj for 
reviews). In the standard concordance ACDM model, 
the overall cosmic fluid contains baryons, cold dark mat- 
ter plus a vacuum energy that fits accurately the current 
observational data and thus it provides an excellent sce- 
nario to describe the observed universe. 

On the other hand, the concordance model suffers 
from, among others [ll|, two fundamental problems: (a) 
The fine tuning problem i.e., the fact that the observed 
value of the vacuum energy density (pA = Ac^/SnG ~ 
10~^'' GeV^) is more than 120 orders of magnitude be- 
low the natural value estimated using quantum field the- 
ory and (b) the coincidence problem 12] i.e., the fact 
that the matter energy density and the vacuum energy 
density are of the same order just prior to the present 
epoch. Such problems have inspired many authors to 
propose alternative candidates to dark energy such as 
A{t) cosmologies, quintessence, A;— essence, vector fields, 
phantoni, tachyons, Chaplygin gas and the list goes on 
(see |13l - [28| and references therein). 

Nowadays, the nature of the dark energy is consid- 
ered one of the most fundamental and difficult problems 
in the interface uniting Astronomy, Particle Physics and 
Cosmology. However, there are other possibilities. For 



instance, one may consider non-standard gravity theo- 
ries where the present accelerating stage of the universe 
is driven only by cold dark matter, that is, with no ap- 
pealing to the existence of dark energy. Such a reduc- 
tion of the so-called dark sector is naturally obtained 
in the so-called f{R) gravity theories 29] (see, however, 
(30|). Even in the framework of the standard general 
relativity theory, is also possible to reduce the dark sec- 
tor by considering the presence of inhomogeneities [sij . 
quartessence models [321 ) . as well as the gravitationally 
induced particle creation mechanism [s^, HI, [38l - l5l| . In 
what follows we focus our attention to the last approach 
by considering the gravitational creation of cold dark 
matter in the expanding Universe. 

The basic microscopic description for gravitational 
particle production in an expanding universe has been 
investigated in the literature by many authors starting 
with Schrodinger [sl]. In the late 1960s, independently 
of Schrodinger's work, this issue was investigated by L. 
Parker and others [SJ, [s^ by considering the Bogoli- 
ubov mode-mixing technique in the context of quantum 
field theory in curved spacetime. The basic physical 
effect is that a classical non-stationary background in- 
fluences bosonic or fermionic quantum fields in such a 
way that their masses become time-dependent (see, e.g., 
36] for more recent works). For a real scalar field in a 
flat Friedniann-Lemaitre-Robertson- Walker (FLRW) ge- 
ometry described in conformal coordinates, for exam- 
ple, the key result is that the effective mass reads [37j , 
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a" /a, where m is the "Minkowskian" 



constant mass of the scalar field, a{ri) is the scale factor, 
and the derivatives are computed with respect to the con- 
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formal time. This time dependent mass accounts for the 
interaction between the scalar field and the geometry of 
the Universe. When the field is quantized, this leads to 
particle creation, with the energy for newly created par- 
ticles being supplied by the classical, time-varying grav- 
itational background. 

Macroscopically, the first self-consistent formulation 
for matter creation was put forward by Prigogine and 
coworkers 13811 and somewhat clarified by Calvao, Lima 
and Waga |39l | (see also [13] ) . It was shown that mat- 
ter creation, at the expenses of the gravitational field is 
also macroscopically described by a negative pressure, 
and, potentially, can accelerate the Universe. Within 
this framework, various interesting features of cosmolo- 
gies with particle creation have been discussed in [4l| - |47| 
(see also [43] for recent studies on this subject). More 
recently, the corresponding effects on the global dynam- 
ics of the particle creation regime has been investigated 
extensively by a number of authors (see j49l - l5ll | and ref- 
erences therein). In particular, it was also found that a 
subclass of such models depends only of one free parame- 
ter and that the evolution of the scale factor is exactly the 
same of ACDM models [5l| . Naturally, in order to know 
if creation of cold dark matter (CCDM) models provide 
a realistic description of the observed Universe, its via- 
bility need to be tested by discussing all the constraints 
imposed from current observations both for background 
and perturbative levels (structure formation). 

In this context, the basic aim of the present work is 
twofold. First, we place constraints on the main param- 
eters of CCDM cosmologies by performing a joint likeli- 
hood analysis involving the shift parameter of the Cos- 
mic Microwave Background (CMB) [3] and the observed 
Baryonic Acoustic Oscillations (BAO) [5], and the lat- 
est SN la data (Constitution) 7]. Secondly, for a wide 
class of matter creation cosmologies, we also develop the 
linear approach for the density perturbations. Analyti- 
cal solutions for the differential equation governing the 
evolution of the growth factor and some properties of 
the large scale structure (cluster formation) are also dis- 
cussed and compared to the predictions of the ACDM 
model. 

The paper is planned as follows. The basic theoretical 
elements of the problem are presented in section 2, where 
we introduce the basic FLRW cosmological equations for 
CCDM cosmologies. In section 3 and 4 we present the 
specific CCDM scenarios and derive the constraints on 
their parameters based on a statistical joint analysis in- 
volving the shift parameter of the Cosmic Microwave 
Background (CMB) the Baryonic Acoustic Oscilla- 
tions (BAO) 0, and the latest SN la data (Constitution) 
. In section 5 we study the evolution of linear perturba- 
tions, while in section 6 we present the corresponding the- 
oretical predictions regarding the formation of the galaxy 
clusters with basis on the Press-Schecther formalism. Fi- 
nally, in section 7 we draw our conclusions. Throughout 
the paper we consider Hq = 70.5 km/sec/Mpc as given 
by the WMAP 5-years data 0]. 



2. CREATION COLD DARK MATTER (CCDM) 
COSMOLOGIES: BASIC EQUATIONS 

The nontrivial cosmological equations for the mixture 
of radiation, baryons and cold dark matter (with cre- 
ation of dark matter particles), and the energy conser- 
vation laws for each component have been investigated 
thoroughly by j40i - t42i . |49- 51] . In this framework, for a 
spatially flat FLRW geometry the basic equations which 
govern the global dynamics of the universe in the matter 
dominated era {prad = 0) are given by 49^511 
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5Hpbar — 0, 
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'iHpdm — Pdm^, (3) 



which implies that 



3Hp,n = pdm^, where p,„ = pdr, 



Pba 



(4) 



In the above set of differential equations, an over-dot 
denotes derivative with respect to time, pbar and pdm, are 
the baryonic and dark matter energy densities, H = a/a 
is the Hubble parameter, whereas Pc, corresponds to the 
creation pressure. The quantity F is the so called creation 
rate of the cold dark matter and it has units of {time)~^. 
The creation pressure is defined in terms of the creation 
rate and other physical quantities. For adiabatic creation 
of cold dark matter, it has been found [H, H^l that the 
creation pressure is given by (see also Prigogine et al. 
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(5) 



Using Eqs. ([TJ and (|4]) we can obtain the following useful 
formula: 



H 



3 H' 



(6) 



Expressions (H])-® show how the matter creation rate, F, 
modifies the basic quantities of Einstein-de Sitter cosmol- 
ogy. Note that if the creation rate is negligible, F << iJ, 
then, the creation pressure is negligible, p,n oc a~^, and 
the solution of the above equation reduces to that of the 
Einstein de-Sitter model, namely, H{t) = 2/3t. 



3. TWO SPECIFIC CCDM MODELS 

Although the precise functional form of T{t) is still 
missing, a number of different phenomenological param- 
eterizations have been proposed in the literature treat- 
ing the time-dependent T{t) function ( [49l - l5l| and refer- 
ences therein). In what follows, we consider two different 
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flat matter creation models namely Lima, Silva & Santos 
(hereafter LSS [i^) and Lima, Jesus & Oliveira (here- 
after LJO [HH) respectively. With the aid of the current 
observational data, we will also put stringent constraints 
on their ft'ee parameters 



3.1. The LSS Model 

In the LSS scenario, the functional form of F is phe- 
nomenologically parameterized by the following linear ex- 
pansion in the Hubble parameter (see LSS for more de- 
tails) 



r = i-yHQ + 3I3H. 



(7) 



It is worth noticing that in its originally formulation by 
LSS, both the baryonic and radiative contributions were 
neglected (pbar = Prad — 0). Such incompleteness was 
cured in a subsequent paper by Steigman et al. [50l |. 
The advantage of the LSS model is that it is analytically 
described and due to this feature we will consider it as an 
interesting and preparatory toy model for more physical 
CCDM scenarios. 

The mass density in LSS reduces to the dark matter 
density, pm = Pdm (see Eq. S]). Without going into the 
details of that model, let us present its main features. 
Based on the latter considerations the basic Friedmann 
equation becomes (see Eq. ^ 



0. 



(8) 



Performing the integration of Eq.® we derive the fol- 
lowing Hubble function: 



H{t) = Ho 



(e^*-l)' 



(9) 



and by integrating the latter we obtain the evolution of 
the scale factor 



a(t) = 



1 



(1 + ^) 



1-7-/3 
7 



3'fHot/2 



1) 



(10) 

Evaluating Ea. ((TO)) at the present time (a = 1 or z = 0), 
we can define the present age of the Universe 
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(11) 



Now utilizing equations © and (ITUl) . we find after 
some algebra, that the normalized Hubble flow takes the 
following form [s^l 



H{a) 
Ho 



7 



1-7-/3 



-3(l-/3)/2 



1-/? 1-/3 



(12) 



We now study the conditions under which an inflection 
point exists in our past, implying an acceleration phase of 



the scale factor. This crucial period in the cosmic history 
corresponds to a{ti) = and tj < to- Differentiating 
twice Eq. (jlOp . we simply have: 



1 



ai 



1 + zi 



27 



(1-3/3) (1-7-/3) 



(13) 



Note that if 7 = 0, there is no transition from an early 
decelerating to a late time accelerating Universe [53|. In 
this case, one can also prove that the expansion of the 
Universe always decelerates for < /3 < 1/3 and always 
accelerates for 1/3 < /3 < 1. The latter condition implies 
that the parameter 7 has to obey the following restriction 
< 7 < 1. When 7 = 0, negative values of /3 means 
destruction of particles, and, for completeness, such a 
possibility will also be considered in the discussion of 
some background tests (see section 4). 



3.2. The LJO Model 

Now, let us consider an alternative CCDM scenario 
that includes the baryonic and radiative components, and 
was also proposed with the intention to solve the cluster 
problem in such a framework [sil ]. In this context, the 
particle creation rate reads 



F = SflA 



Pec 
Pdn 



H, 



(14) 



where Q\ (called a in the 51]) is a constant, pco = 
3Hq I SttG is the present day value of the critical density, 
and the factor 3 has been maintained for mathematical 
convenience. We stress that such CCDM scenario mimics 
the global dynamics of the traditional A cosmology fsii] . 
Indeed, the combination of equations ([Ij, ([6]), and 
leads to the following formula: 



h + Ih' = 1^ahI 



a solution of which is given by 
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Using now the deflnition of the Hubble parameter H = 
d/a, the scale factor of the universe a(t), normalized to 
unity at the present epoch, evolves with time as: 
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(17) 



where f2m = 1 — f^A = 1 — a, can be viewed as the effective 
matter density parameter at the present time. For an 
arbitrary redshift, Clrn.iz), measures the mass density that 
is effectively clustering [5l|. The cosmic time is related 
with the scale factor as 



tia) 



-sinh 




(18) 
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Combining the above equations we can define the normal- 
ized Hubble expansion as a function of the scale factor: 



Eia) 



Hja) 
Ho 



-3 



1/2 



(19) 



Obviously, the LJO model contains only one free param- 
eter. The inflection point [namely, the point where the 
Hubble expansion changes from the decelerating to the 
accelerating regime, a{ti) — 0] takes place at: 



ti = 



-sinh 




ai 



20^ 



-I 1/3 



(20) 



In what follows we constrain the free parameters of such 
models coming from the background tests. For the sake 
of clarity, the study of the perturbed models and the asso- 
ciated constraints coming from galaxy cluster formation 
will be postpone to sections 5 and 6. 



4. BACKGROUND TESTS: BAO, CMB SHIFT 
PARAMETERS AND SNE lA DATA 

In this section we briefly present the statistical analy- 
sis used in order to constrain the matter creation models 
presented in the previous section. First of all, we consider 
the Baryonic Acoustic Oscillations (BAO). This kind of 
estimator is produced by pressure (acoustic) waves in the 
photon-baryon plasma in the early universe, generated 
by dark matter overdensities. Evidence of this excess 
was recently found in the clustering properties of the lu- 
minous SDSS red-galaxies 5j and it provides a "stan- 
dard ruler" with which we can constrain dark energy 
models. In particular, we use the following estimator: 

2/3 



J as a 



da 
ma) 



measured from the 



[z2£;(a,)]i/3 

SDSS data to be A = 0.469 ± 0.017, where 2, = 0.35 [or 

Us = {I + Zs)~^ ~ 0.75]. Therefore, the corresponding 

,2 



Xbao function is simply written 



Xbao(p) 



[A{p) - 0.469]^ 
0.0172 



(21) 



where p is a vector containing the cosmological parame- 
ters that we want to fit. 

On the other hand, a very accurate and deep geomet- 
rical probe of dark energy is the angular scale of the 
sound horizon at the last scattering surface as encoded 
in the location if^ of the first peak of the Cosmic Mi- 
crowave Background (CMB) temperature perturbation 
spectrum. This probe is described by the CMB shift pa- 
rameter [H, [U and is defined as i? = V^m J^^ a^'E(a) ■ 
The shift parameter measured from the WMAP 5-years 
data 0] is i? = 1.71 ± 0.019 at zis = 1090 [or ais = 
(1 -I- zis)'^ ~ 9.17 X 10^"*]. In this case, the X^mb func- 
tion is given 



Xcmb(P) 



[Rjp) - 1-71]^ 

0.0192 



(22) 



Note that the measured CMB shift parameter is some- 
what model dependent but mostly to models which are 
not included in our analysis. For example, in the case 
where massive neutrinos are included or when there is 
a strongly varying equation of state parameter. The ro- 
bustness of the shift parameter was tested and discussed 
in [54j]. 

Finally, we use the 'Constitution set' of 397 type la 
supernovae of Hicken et al. [7]. In order to avoid possi- 
ble problems related with the local bulk flow, we use a 
subsample of the overall sample in which we select those 
SNIa with z > 0.023. This subsample contains 351 en- 
tries. The corresponding xInj^ function becomes: 



XSNIa(p) 



351 

E 



M*^ai,p)-/x°''^(a,) 



(23) 



where ai = (1 -I- z^)^^ is the observed scale factor of the 
Universe, Zi is the observed redshift, /i is the distance 
modulus fi — m — M — 5\ogdL + 25 and dL(a,p) is 
the luminosity distance dL{a,p) = ^ /J y2H{y) where c 
is the speed of light. We can combine the above probes 
by using a joint likelihood analysis: £tot(p) = -Cbao x 

/^cmb X £sNIa Or X?ot(P) = XbAO + Xcmb + XSNIa m O^der 

to put even further constraints on the parameter space 
used"^. 

Performing now our statistical analysis we attempt 
to put constrains on the free parameters. In partic- 
ular, we sample /3 £ [-0.2,0.2] and 7 G [0.1,1] in 
steps of 0.001. Thus, the overall likelihood function 
peaks at ;9 = -0.15 ± 0.007 and 7 = 0.86 ± 0.01 with 
XtotiPjl) — 432.5 (rfo/=351). Using the latter cosmo- 
logical parameters the corresponding current age of the 
universe is found to be to ^ 14.8 Gyr while the inflec- 
tion point is located at ai ~ 0.44. This corresponds to 
a redshift zj ~ 1.26, which is substantially higher than 
in the case of the concordance model. If we now impose 
fH — 0, the SNIa likelihood analysis provides a best fit 
value of 7 = 0.66 ± 0.02, in agreement with [5^. Note 
that [H^] used only the SNIa data. We find that although 
our SNIa likelihood analysis provides a similar 7 value to 
that found in fEd], our joint likelihood funtion peaks at 
7 = 0.41 ± 0.01 but with a poor quality fit: the reduced 
xiot is =i 3. 

This simply means that the functional form E{a) = 
7 -I- (1 — 7)a~'^/2 is unable to fit the observational data 
simultaneously at low and high redshifts. We confirm 
this point by using the CMB shift parameter only, finding 
that the corresponding likelihood function peaks at 7 ~ 
0.20. This value is ~ 3.4 times larger than that provided 



^ Likelihoods are normalized to their maximum values. We will 
report Icr uncertainties of the fitted parameters. The overall 
number of data points used is Ntot = 353 and the degrees of 
freedom: dof= Ntot — nat , with ngt the number of fitted param- 
eters, which varies for the different models. 
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by the SNIa+BAO solution 7 ~ 0.68 (SNIa only points 
7 ~ 0.66). This fact alone suggests that for the class 
of models with a constant creation rate F — S-fHo, are 
unable to provide a quality fit of the basic cosmological 
data in all the relevant redshift ranges. 

Now, let us constrain the LJO model with the ob- 
servational data from background tests. From our joint 
analysis involving BAO + CMB Shift Parameter + SNe 
la data, we find that the best fit value (within la un- 
certainty) is nm = 0.28 ± 0.01 with Xtoti^m) ^ 431.5 
((io/=352). This result is in good agreement with re- 
cent studies [iHl, 1^, 0] of the concordance A cosmology. 
Therefore, the current age of the universe is to ~ 13.9Gyr 
while the inflection point is located at a/ ~ 0.58 (hence 
at redshift zj ~ 0.72). 

Summing up, although solving the acceleration and age 
problems, the LSS type scenarios is endowed with severe 
difficulties even for the set of background tests analysed 
here. On the other hand, the LJO scenario provides good 
fits for SNe + BAO + CMB shift parameter. A more 
detailed analysis on this issue (background tests) it will 
be published in a forthcoming paper. 



5. MATTER DENSITY PERTURBATIONS IN 
CCDM COSMOLOGIES 

Let us now derive the basic equation that governs the 
evolution of the mass density contrast, modeled here as 
an ideal fiuid in a 'quasi'-Newtonian or neo-Newtonian 
framework [ll]. In this approach, background equations 
are formulated in a way that the isotropic pressure be- 
comes dynamically relevant even for FLRW cosmologies 
where all the physical properties may depend only of 
time. This allows us to describe an accelerated expansion 
of the Universe powered by an effective negative pres- 
sure in a Newtonian framework. Naturally, even consid- 
ering that the neo-Newtonian approach reproduces the 
GR background dynamics exactly, small differences may 
occur at the perturbative level. However, as discussed 
by Reis [llj, the GR first-order perturbation dynamics 
and its neo-Newtonian counterpart coincide exactly for 
a vanishing sound speed. Indeed, for constant equations 
of state it has been demonstrated that the correct large- 
scale behavior in the synchronous gauge are also repro- 
duced. This should also be expected since the observa- 
tional data correspond to modes that are well inside the 
Hubble radius, for which the use of a neo-Newtonian ap- 
proach seems to be adequate (for recent papers in this 
subject see [S^l and references therein) 



If) +V..(pu)=pF, 



du 

'dt 



(u • V,.)u 



and 



V^?$ = 47rGp-Ae//, 



(24) 



(25) 



(26) 



where (r, t) are the proper coordinates, u is the velocity 
of a fiuid element of volume, p is the mass density and 
$ is the gravitational potential. Since we are working 
within the context of cosmological models with matter 
creation we have modified the continuity equation ([M)) 
by including the "standard" source term. Utilizing Eq. 
^ and Eq. ([5]) , the quantity Aef / is defined as 



A 



A - UttGpc a + ^nGpdm 



F 



(27) 



where for completeness and future comparison to the 
present cosmic concordance model ( ACDM) , we have also 
included the cosmological A-term. 

Now, by changing from proper (r, t) to comoving (x, t) 
variables, the fiuid velocity becomes |58| 



u = ax + ax = ax + v(x. t), 



(28) 



while the corresponding differential operators take the 
following forms 



S/x = V = aV,. 



and 



d_ 

dt 



(29) 



(30) 



where x = r/a and v(x, i) is the peculiar velocity with 
respect to the general expansion. Note that the mass 
density is written as 



Pm{t)[l+S{ji,t)]. 



(31) 



In this context, using eqs.(|l]), ([^ and (PU)) and neglect- 
ing second order terms (5 <C 1 and u <C it) it is routine 
to rewrite eqs.dMl), (ES]) and ^ 



dv 

ax + — + 
ot 



1 

a 



(32) 



5.1. Basic Formalism 

In what follows we will adopt a nonrelativistic descrip- 
tion with basis in an extended continuity equation to- 
gether the Euler and Poisson equations. In virtue of the 
particle creation process they take the form: 



Vv 



dt pm 



-^V2$ = 47rGp„(l + <5)- Ae//. 



(33) 



(34) 
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Following the notations of [58| , we write the gravitational 
potential as 

2 1 

$ = (/)(x, t) + -nGp,na^x^ - -Aeffa^x^. (35) 

o D 

Thus utilizing Eqs.(IT|), ^ and (|27| we derive after some 
algebra that 



and 



5v V0 

7^+^"^ = ' 

at a 



(36) 



(37) 



Finally, by taking the divergence of Eq. (j36|) and using 
Eqs. p3p and p7p . we obtain the time evolution equation 
for the growth factor D{t) = S/A{x.) [where A(x) is an 
arbitrary function] 

D + {2H + Q)b - (47rGp,„ - 2HQ - D = 0, (38) 

where 



Qit) = 



Pm 



(39) 



Now, it is clear how the matter creation term affects the 
growth factor via the function Q{t). 

It should be noticed that the background creation pres- 
sure affects both the global dynamics as well as the fluc- 
tuations via Q{t). In the above results, it was also implic- 
itly assumed that the produced particles have negligible 
velocities as measured by the observers in the comov- 
ing framework which implies that any additional internal 
pressure term is negligible. In fact this is a reasonable 
assumption since we are working with a non-relativistic 
phase of universe expansion (matter era). The second 
assumption is that the matter production is strictly ho- 
mogeneous which implies that possible contributions of 
Q{t) at the perturbative level is practically zero. Thus, 
within the context of the latter assumptions, the advan- 
tage of employing a neo-Newtonian approach is a gain in 
simplicity and transparency of all computations. 

In order to solve the above differential equation we 
need to know the functional form of the quantity Q{t) (or 
r). In the case of a negligible matter creation rate, F << 
H [Q{t) ~ 0], the above equation reduces to the usual 
time evolution equation for the mass density contrast [11] 
for which the growth factor is DEds{o,) = o,- As one may 
check, by solving Eq. ([55)1 for the A cosmology [Q{t) — 0] , 
we also derive the well known perturbation growth factor 
(see H): 



DAia) = 



x'iE^ixY 



(40) 



where Q,m = 1 — ^^a is the matter density parameter and 
E{a) = H{a)/Ho is the normalized Hubble function 



E{a) = {ilA + n^a-^) 



(41) 



It is also interesting to mention here that one can ex- 
plicitly derive Eq. p8|) in the framework of cosmological 
models with a time varying vacuum energy density [59|. 
Note, that for simplicity throughout the rest of the pa- 
per we will use geometrical units (c = SttG =1). In the 
analytical treatment below, unless explicitly stated, we 
consider A = 0. 



5.2. Growth of Fluctuations in CCDM Models: 
Analytical Solutions 

Let us now discuss thoroughly the time evolution equa- 
tion of the mass density contrast in the linear regime in 
order to test the implications of CCDM models on the 
structure formation process. In order to illustrate our re- 
sults with a basic application, we also consider later the 
formation of galaxy clusters through the so-called Press- 
Schetcher formalism. 



5.2.1. Fluctuations in LSS Model 

Using the time evolution equation for the mass den- 
sity contrast as given by ([38]), we now derive the growth 
factor of fluctuations for LSS cosmology. To begin with, 
we change variables from t to a new one defined by the 
following transformation 



y 



= e"*/^ - 1 with n = 3-/Ho. 



(42) 



In this context, using equations (P), (g]), dH), (fTO|) . 
and (H^ we obtain 



H{y) 



Pmiv) 



Q{y) 



Qiv) 



njy + l) 
My 

9 9 ' 

Mr 

^ njSy - (m - 3)] 

My 

n^ip~3)iy+l) 
2py^ 



(43) 



where p — 3{1 — /3). We can now rewrite Eq. ((38)) as: 

p^iy + l)y^D + pyf{y)D + 2g{y)D = 0, (44) 
where prime denotes derivatives with respect to y and 



f{y) - {lQ + p){y + l)^ip 
g{y) = 9{y + 1) - 7p + p' . 



(45) 



Factors of n drop at the end of the calculation. Notice, 
that this variable is related with the scale factor as 



y-^a^^^. 

M - 37 



(46) 
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We find that eq. p4|) has a growing solution of the form 



where 



3 \/7 \/A , VA 

1^1 = + ^r~ ^2^1 + — , 

A ==28 -4^ + //. 



(47) 



(48) 



Notice that the quantity F is the hypergeometric function 
and C is a constant. Inserting ea. (P5| into cq. (|47p . we 
obtain the growth factor D{a) as a function of the scale 
factor: 

(49) 



where 



37 



^ - 37 



(-10+3ai+V^)/2m 



(50) 



5.2.2. Fluctuations in LJO Models 

Let us now proceed further to solve analytically the 
differential Eq. (pS)) in order to investigate the matter 
fluctuation field of the LJO model in the linear regime. 
To do so, we change variables from f to y according to 
the transformation 



y — coth t 



(51) 



In this framework, the differential equation (|52p can 
be brought into the following expression 



dD 



6x^{x + iy -r^+xq{x) —~b{x)D :^0, (55) 



dx 



q{x) = (x + l)(5a; - 4) b{x) = {x^ - bx - 3), (56) 
where 



x = y^-l= {^\ a 



x+1 = y^ 



E\a) 



, (57) 



in which we have used equations (fT6|) , (|17l) and (IT9t . 
The solution of the main differential equation ([55)) is 



D{x) =C- 



,(5-77)76 
X + 1 



5 7, 

+ 7:,2(k+ -),-a; 
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, (58) 



where k = — (4+-\/7)/6 andC is a constant. Inserting Eq. 
(1571) into Eq. (|58p . we finally obtain the growth factor 
_D(a) as a function of the scale factor: 



D{a) ^ Ci 



where 



-(5-v7)/2 



K, K H — ,2(kH — , — - — a 
6' ^ 6^' 

(59) 



(5-\/7)/6 



(60) 



We would like to stress here that the above solution of 
the growth factor is valid when 17 a 7^ (or i7„i 7^ !)• 



Using ([5T|) we find, after some calculations, that equation 
takes on the form 



3y2(j;2 _ 1)2^" + 2yj{y)Fl - 2g{y)D = 0, (52) 
where 

fiv) - - i)(y' - 3), 5(y) = (y'^ - V + 3). (53) 

Note that factors of Hq drop at the end of the calculation. 

In deriving Eq. (j52p we have substituted the various 
terms in ([55| as a function of the new variable [see Eq. 
([5T|) ]. For instance, from equations dU, (flB)) and ([M|) . we 
have 



i7(y) = V^A^oy, 
370^^0 



Q(y) 



y 



Q(y) 



9VriAgo (y' - 1) 

2j/2 



(54) 



5.3. Linear Growth: Analysis and Summary 

Let us now discuss the evolution of the growth factor 
as a function of the redshift [D{z), z = a^^ — 1] as pre- 
dicted by the two CCDM cosmologies (see Eqs. (HH]) and 
([551) ) . For comparison we also present the growth factor 
provided by the standard ACDM cosmology (see Eq. l40l) . 

In figure 1 we display the basic results. First, we no- 
tice that the growth factors are normalized to unity at 
the present time [z = 0). It is also evident that the 
LSS growth factor (dashed line) is much greater with re- 
spect to the other 2 models, LJO (dot line) and A (solid 
line). In this framework, we find that for z < 0.8 the 
LSS growth factor (dashed line hereafter LSSi) starts to 
decay, implying that cosmic structures cannot be formed 
via gravitational instability. The same general behavior 
seems to hold also for the LSS model with (3 — Q and 
7 = 0.66 (dot-dashed line-hereafter LSS2). 

Let us now compare the predictions of LJO and ACDM 
growth factors. First we stress a remarkable result: de- 
spite the fact that the two models share the same global 
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FIG. 1: The evolution of the growth factor. The hnes 
correspond to LSSi [dashed, (/?,7) = (-0.15,0.86)], LSS2 
[dot-dashed, (/3,7) = (0,0.66)], LJO (dot) and A cosmology 
(solid). Note that the growth factors of all models were scaled 
to unity at the present time. Despite the fact that LJO and 
ACDM models share the same global dynamics, wee see that 
they trace differently the evolution of the matter fluctuation 
field (compare the dotted and solid lines). 



dynamics (compare Eqs. and (HT|) ). they trace dif- 
ferently the evolution of the matter fluctuation field. In 
particular, close to the present epoch (z < 0.4) the LJO 
growth factor (dot-dashed line) reaches almost a plateau, 
which means that the matter fluctuations are effectively 
frozen. Between 0.4 < z < 1.6, the growth factor in 
the LJO model is greater than that of the concordance 
ACDM cosmology. In particular, the corresponding de- 
viation {1-D/Da)% lies in the interval [-35, -19]%. As 
an example, prior to the inflection point {zj ^ 0.72) we 
find - -28%. 

It is worth noting that in the interval of 1.6 < z < 2.8 
the LJO growth factor tends to the A solution. In- 
deed, at the epoch of z ^ 1.62 (a ~ 0.38), in which the 
most distant cluster has been found [60| . the deviation 
(1 - D/Da)%, of the growth factor D{z) for the LJO 
scenario with respect to the A solution D\{z) is —5%. 
To this end, for z > 2.8 we find that D{z) < Da{z) and 
the corresponding deviation {1 — D/Da)%, remains close 
to ~ 28%. Naturally, one may expected that the differ- 
ences among the growth factors will affect the predictions 
related with the formation of the cosmic structures (see 
next section). In particular, such results point to the pos- 
sibility of a crucial observational test confronting CCDM 
with ACDM scenarios in a perturbative level. 



THE FORMATION AND EVOLUTION OF 
COLLAPSED STRUCTURES 



In this section we study the cluster formation pro- 
cesses in CCDM cosmologies by using the standard Press- 
Schecther formalism (6l| . In such an approach, the be- 
havior of the matter perturbations is described by as- 



suming that the density contrast behaves like a Gaussian 
random field. As it is widely known, the cluster distribu- 
tion basically traces scales that have not yet undergone 
the non-linear phase of gravitationally clustering, thus 
simplifying their connections to the initial conditions of 
cosmic structure formation. 

The basic aim here is to estimate the fractional rate of 
cluster formation (see 62, 63]). In particular, these stud- 
ies introduce a methodology which computes the rate at 
which mass joins virialized halos (such as galaxy clusters) 
which grow from small initial perturbations in the uni- 
verse, with matter fiuctuations, 6, greater than a critical 
value 6c- 

Now, by assuming that the density fiuctuation field, 
smoothed at the scale R (corresponding to a mass scale of 
M = AirpR^ /3, with p the mean background mass density 
of the Universe) , is normally distributed with zero mean, 
then the probability that the field will have a value 6 at 
any given point in space is: 



r{s,z) 



1 



2'Ka{R,z 



■exp 



2a2(i?, 



(61) 



where the variance of the Gaussian field, a^{R,z), is 
given by: 

1 f°° 

a^{R,z)^— k^P{k,z)W^{kR)dk, (62) 
27r2 Jo 

with P{k, z) the power spectrum of density fluctuations 
which evolves according to P{k, z) = P{k, 0)D^{z), with 
D{z) the growing mode of the density fiuctuations evo- 
lution, normalized such that D{0) = 1. 

Finally, W{kR) is the top-hat smoothing kernel, given 
in Fourier-space by: 



W{kR) 



[kpy 



[sm{kR) - kRcos{kR)] 



(63) 



We can now estimate what fraction of the Universe, at 
some reference redshift z and above some mass thresh- 
old M, has collapsed to form bound structures. To this 
end we need to integrate the probability function given 
by Eq. ([5T|) over all regions that at some prior redshift 
had overdensities which by the reference redshift have 
increased to above the critical value, Sc{z), which in an 
Einstein-deSitter universe is ~ 1.686 and varies slightly 
for different values of r2,„ f63,'6^. Therefore, this fraction 
is given by [63],[66],^67.],,68j: 



F{M,z)^ / V{5,z)d5, 

•^*c(z) 



(64) 



and performing the above integration, parametrizing the 
rms mass ffuctuation amplitude at i? = 8 Mpc which 
can be expressed as a function of redshift as ^(M, z) — 
aii{z) — D{z)as, we obtain: 



1 - erf 



\/2a8(z) 



(65) 
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FIG. 2: The predicted fractional rate of cluster formation as 
a function of redshift for erg — 0.80. The lines represent: (a) 
LSSi [dashed, (/3,7) = (-0.15,0.86)], (b) LSS2 [dot-dashed, 
(^,7) = (0,0.66)], (c) LJO (dot) and (d) A cosmology (solid). 
Note that the clusters in LSS type models form earlier {z ~ 
7 — 10) in comparison to those produced in the framework of 
the LJO and ACDM models (2 ~ 2). 

It is worth noticing that the above generic form of Eq. 
is heavily dependent on the choice of the background 
cosmology and the matter power spectrum normaliza- 
tion, as ■ For the sake of comparison, in what follows we 
consider two values erg. The first is the WMAP5 result 
CTg — 0.80 m whereas the second is erg — 0.95. 

Finally, we normalize the above probability to give the 
fraction of mass in bound structures which have already 
collapsed by the epoch z, divided by the corresponding 
fraction in structures which have collapsed at the present 
epoch (z = 0), 

F{z) ^ T{z)/T{0) . (66) 

In Fig. 2, by assuming the matter power spectrum nor- 
malization, (jg — 0.80, we present in a logarithmic scale 
the behavior of normalized structure formation rate as a 
function of redshift for the present cosmological models. 

In the context of the LSSi matter creation model (see 
dashed line), we find that prior to z ~ 0.8 (or z ~ 1.1 for 
the LSS2 model, dot-dashed line) the cluster formation 
has effectively terminated due to the fact that the matter 
fluctuation field, D{z), effectively overs (see section 5.3). 
Also, the large amplitude of the LSS fluctuation field 
(see also Fig. 1) implies that in this model galaxy clus- 
ters appear to form much earlier (z ^ 7) with respect 
to the LJO (dot-dashed) and A (solid line) cosmologi- 
cal models. Indeed, for the latter cosmological models 
we find that galaxy clusters formed typically at z ^ 2.2. 
From the observational point of view, it is interesting to 
point here that the most distant cluster has been found 
at z ~ 1.62 [iOj. We would like to stress that the LSS2 
predicts an early rapid cluster formation that takes place 
at the redshift of reionization (z ~ Zreion ~ 10), at which 
the universe was reionized from the neutral state to a 
fully ionized state. Obviously, this result rules out the 




z 



FIG. 3: The predicted fractional rate of cluster formation as 
a function of redshift for using a higher value for the nor- 
malization of the power spectrum, as — 0.95. The meaning 
of the various lines is the same of Fig. 2. We see that for a 
larger value of ug the corresponding cluster formation rate is 
displace to higher redshifts (compare to Fig. 2). 

LSS2 model. 

We also stress that the amount of cluster abundances 
between the two models remains almost the same for 
1.6 < z < 2.2. This is to be expected due to the fact 
that the two cosmological models have almost the same 
growth factors (see also section 5.3). Then for z < 1.6 
the LJO cluster formation rate becomes larger than that 
of the concordance model and finally it terminates for 
z < 0.4. 

In Fig. 3 we also display the normalized structure for- 
mation rate as a function of redshift, but now for a higher 
value of the power spectrum parameter (erg = 0.95). 
Comparing Figs. 2 and 3, we see that the basic effect 
of a larger value of ag is that the corresponding cluster 
formation rate moves to higher redshifts. 



7. CONCLUSIONS 

In this work, we have investigated (analytically and 
numerically) the overall dynamics of two cosmological 
models in which the dark matter creation process pro- 
vides the late time accelerating phase of the cosmic ex- 
pansion without the need of dark energy. Such scenarios 
termed here by LSS [4^] and LJO ^l[ are phenomenolog- 
ically characterized by two distinct creation rates, F (see 
section 3). In the first scenario (LSS), the creation rate is 
r = 3jHo+3/3H while in the second one it is proportional 
to the Hubble parameter, namely: F = SilAiPco/ Pdm)H , 
where pco = SHq/SttG is the present day value of the 
critical density. 

It should also be stressed that the phenomenological 
approach adopted here cannot determine the mass of the 
particles, as well as whether their nature is fermionic or 
bosonic. In order to access the mass of the particles, and, 
therefore, the nature of the CDM particles, it is necessary 
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to determine the creation rate, F, from quantum field the- 
ory in FLRW space time. In principle, such a treatment 
must somewhat incorporate a source of entropy since the 
matter creation process is truly an irreversible process. 
In particular, in the case of adiabatic matter creation 
considered here, both the entropy (S) and the number 
of particles (N) in a comoving volume increase but the 
specific entropy (S /N) remains constant 39, 41^, ■ 

In this context, by using current observational data 
(BAOs, CMB shift parameter and SNIa) we have first 
performed a joint likelihood analysis in order to put tight 
constraints on the main cosmological parameters. Subse- 
quently, trough a linear analysis we have also studied the 
growth factor of density perturbations for both classes of 
creation cold dark matter models, and, finally, by using 
the Press-Schecther formalism we have discussed how the 
cluster formation rates evolve in such scenarios. 

The result of our joint statistical analysis shown that 
the fit provided by the LSS model with 7 — 0.66 and 
(3 = turns out to be of poor quality because it is unable 
to adjust simultaneously the observational data at low 
and high redshift. This result confirm the conclusion by 
Steigman et al [H^ using only SNe la and the determina- 
tion of Zeq, the redshift of the epoch of matter radiation 
equality. For the LSS scenario, we also find that the am- 
plitude and the shape of the linear density contrast are 
significantly different with respect to those predicted by 
the LJO and A models (see Fig. 1). In particular, for 
z < 0.8 the matter fiuctuation field of the LSS model 
practically decays thereby implying that the correspond- 
ing cosmic structures cannot be formed via gravitational 
instability (see also Figs. 2 and 3). 

On the other hand, the creation cold dark matter sce- 
nario (termed here as LJO) provides good quality fits 
of the cosmological parameters at all redshifts and it re- 



sembles the global dynamics of the concordance ACDM 
cosmology by including only one free parameter. Despite, 
the latter the corresponding growth factor has the follow- 
ing evolution with respect to that of the usual A cosmol- 
ogy: (i) prior to the present epoch (z < 0.4) the evolution 
of the LJO growth factor reaches almost a plateau, which 
means that the matter fiuctuations are effectively frozen, 
(ii) between 0.4 < z < 1.6, the growth factor in the LJO 
model is greater than that of the concordance A cosmol- 
ogy, while for 1.6 < z < 2.8 the two growth factors have 
converged and (iii) for z > 2.8 we find that D{z) < Da{z) 
and the corresponding deviation (1 — D/ D\)%, remains 
close to ~ 28%. 

Summarizing, in the case of LSS scenario, the large 
scale structures (such as galaxy clusters) form earlier 
(z ~ 7 — 10) with respect to those produced in the frame- 
work of the LJO and ACDM models (z ^ 2). Therefore, 
in view of the recent cluster observational data, the lat- 
ter are much more favored as compared to the former. 
Our basic conclusion is that the LJO creation cold dark 
matter model which has only one free parameter passes 
all the tests considered here. However, new constraints 
from complementary observations need to be investigated 
in order to see whether this kind of scenario provides a 
realistic description of the observed Universe. 
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